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The first penetration field Hj(T) of Josephson vortices is measured through the onset of microwave 
absorption in the locked state, in slightly overdoped Bi2Sr2CaCu20g+i single crystals (T c « 84 K). 
The magnitude of Hj(T) is too large to be accounted for by the first thermodynamic critical field 
Hd(T). We discuss the possibility of a Bean- Livingston barrier, also supported by irreversible 
behavior upon flux exit, and the role of defects, which relates Hj(T) to the c-axis penetration depth 
A C (T). The temperature dependence of the latter, determined by a cavity perturbation technique 
and a theoretical estimate of the defect-limited penetration field are used to deduce from Hj(T) the 
absolute value of A c (0) = (35 ± 15) jim. 

PACS numbers: 47.32Cc, 74.72h, 78.70Gq 



I. INTRODUCTION 

The phenomenological Lawrence-Doniach model is 
generally used to describe a stack of Josephson-coupled 
superconducting layerstrB. This interlayer Josephson 
tunneling has been established experimentally by dc or 
ac Josephson cfEcct experiments in numerous high-T c 
superconductorsu and is proposed as a candidate mecha- 
nism for superconductivity!! Such discrete layered struc- 
tures have some striking incidence on many properties: 

i) Josephson vortices appear for field parallel to the 
layers, and in case their penetration in_this quasi-2D sys- 
tem is impeded by a surface barrieiQ, the penetration 
field, henceforward noted H^ D (T), is simply inversely 
proportional to the c-axis penetration depth X c (Tp, un- 
like isotropic superconductors (where it is of the order 
of the thermodynamic critical field). The occurence of 
such a barrier was discussed mostly in the framework 
of low-field magnetization measurements performed in 
fields parallel to the layers in NdCeCuOB, Tl-22010 and 
Bi-2212E]. The quantitative estimates p£ A C (T) deduced 
from these data were however disputedliil. 

ii) A C (T) is directly related to the critical current den- 
sity between the layers, Jo(T), and is inversely-propor- 
tional to the Josephson plasma frequency w ps E3tl Both 
quantities ough t to b e discussed within the same theoret- 
ical backgrounduJllj. Direct determination of the plasma 
frequency was performed through infrared reflectivity 
measurements iii_La2_xSr x Cu04li3, YBa2Cu306+<5", Tl- 
2212 or Tl-220lE£l and from microwave absorption mea- 
surements in underdoped Bi-2212 and Bi-2201c3. A large 



body of literature reported a sharp microwave absorption 
line in presence of a static field applied parallel to the 
c-axisEIfEi This absorption line was assigned to Joseph- 



son plasma resonance, whose frequency is r 
the field-dependent interlayer phase coherence 



ed by 
. How- 



ever, this interpretation is still controversialc-l Although 
the geometry of the experiments reported here is differ- 
ent (the external field is parallel to the ab planes), the 
specific field dependence of A c or u> ps may be involved, as 
discussed elsewhereEjlia. 

Therefore, an independent measurement of the abso- 
lute value of A c (in zero applied field) is of interest. To 
date, all the above mentioned properties have been stud- 
ied separately. It is the aim of this paper to bring to- 
gether two different microwave measurements in order to 
obtain the absolute value of A C (T): i) the first penetration 
field of Josephson vortices is measured and shown tp-.be 
related to A C (T), (ii) a cavity perturbation techniqueE3 is 
used to detepnine the temperature variation of AA a f,(T) 
and AA C (T)H. 

In the present paper, we focus mainly on the investiga- 
tion of the penetration of Josephson vortices through sur- 
face resistance measurements at high frequency (10 GHz) 
in Bi-2212c3. The onset of microwave absorption allows 
to determine the penetration field Hj(T) of Josephson 
vortices at different temperatures. The magnitude of 
Hj(T) and the irreversible behavior of the dissipation 
with respect to flux entry and flux exit point at first 
sight toward a Bean-Livingston surface barrier imped- 
ing the penetration of Josephson vortices. However, a 
closer quantitative investigation, which includes the ex- 
perimental determination of the variation AA C (T) of the 
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c-axis penetration depth, and the theoretical calculation 
of the penetration field in the presence of edge or surface 
defects, leads us to the conclusion that Hj(T) is eventu- 
ally controlled by such surface irregularities. Relying on 
these theoretical estimates, we deduce from Hj(T) the 
absolute value of the c-axis penetration depth. 



II. EXPERIMENT 

Microwave dissipation measurements were performed 
in various (generally slightly over-doped) BSCCO sin- 
gle crystals shaped into rectangular platelets of approx- 
imate size a x 6 x c ~ 2 x 1 x 0.03 mm 3 : sample A, 
T c = 86 K, has a transition width AT C « 3 K (as de- 
termined from the range over which the microwave ab- 
sorption drops from normal to superconducting state val- 
ues), sample B with T c = 84 K, AT C m 3 K, sample C, 
T c = 89 K, AT C « 3 K. Two other similar samples (D 
and E) were used for checking the onset of microwave dis- 
sipation with respect to the surface quality as discussed 
below. Finally, the temperature dependence of the pen- 
etration depth was measured in a set of similar samples 
by a cavity perturbation technique at 10 GHz and ac- 
susceptibility at 100 kHz. The details of these measure- 
ments will be discussed elsewherecj while here we shall 
only make use of the temperature variations of AX a b(T) 
and AA C (T). An example of the temperature dependence 
of the surface resistance R s (T) in the a6-planc of slightly 
overdoped (T c = 84 K) BSCCO single crystal is shown 
in Fig. 1 . The extrapolation of this curve to T = (inset 
of Fig. 1) yields estimate R ICS w 110 [ifl, which is, to the 
best of our knowledge, the lowest value ever obtained in 
BSCCO single crystals at 10 GHz. The inset of Fig. 1 dis- 
plays also the linear change with temperature (T < 50 K) 
of AA afc (T) = X a b{T) - X ab (5K). This linear variational 
low T was previously observed in optimally dopedt§~E£l 
and slightly overdopedE2l BSCCO single crystals. Both 
of the above mentioned parameters of the sample sug- 
gest that the quality of the cuprate planes is fairly high. 
We note that the slope of AA a &(T) in the inset of Fig. 1 
is fairhzJaxge (w 25 A/K). It could be a consequent 
dopingcj'u with respect to optimally doped crystalsc 

All samples from different batches exhibit very similar 
properties as far as the magnitude and temperature de- 
pendence of the field penetration is concerned, so that 
experimental results are only displayed for sample A. 
AA C (T) differs among samples with different T c . We shall 
only make use of the data on the samples with the same 
T c (±1 K). 

The experimental set-up was described elsewhereEa. It 
is used to measure the microwave losses as a function of 
the applied magnetic field (0-100 Oe) and temperature 
(50-90 K, measurements at temperatures lower than 45 K 
are hindered by the increasing noise of the set-up) . 

The microwave field hi lies within the afo-plane, so that 
the induced microwave currents flow both within the ab- 
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FIG. 1. Surface resistance R S (T) in the a6-plane of slightly 
overdoped BSCCO single crystal. The inset shows the low 
temperature behavior of Rs{T) and of the penetration depth 
AA a6 (T). 



plane and along the crystallographic c-axis. The static 
magnetic field H is applied in the afe-plane perpendicular 
to the microwave field. A computer-controlled goniome- 
ter allows to select its orientation 8 with respect to the 
afr-plane. To locate the 8 — position, we take advantage 
of the lock-in transition evidenced earliero. The set-up 
measures the variation of the power dissipated in the cav- 
ity as the magnetic field is swept at fixed temperature, 
hence yields the field inducecLimaginary part x"(H) of 
the macroscopic susceptibility^ (as long as the dissipa- 
tion is ohmic, the so-called linear regime). This latter 
point has been checked for all the data shown hencefor- 
ward. 



III. RESULTS 

Figure 2 displays the change of dissipation x"(H) — 
x"(0), starting from zero field (within ±0.1 Oe) mea- 
sured in sample A for various orientations of the ap- 
plied field close to the a6-plane: 0° < 8 < 3° (only 
the 0° and 2° are displayed in Fig. 2) and in a low-field 
range: < H < 25 Oe, at three typical temperatures 
(T — 78 K, 65 K, 50 K). 

After each field sweep, the sample was warmed through 
T c and then cooled again in zero field, in order to avoid 
any possible vortex pinning when studying the penetra- 
tion starting from zero field. The dissipation of Joseph- 
son vortices is characterized by the fact that it does not 
depend on the angle (Fig. 2), as long as these vortices re- 
main locked. According to our previous study, the dissi- 
pation regime displayed in Fig. 2 comes only from locked 
Josephson vorticesEj and holds up to ~ 30 Oe. 

As the field increases, an onset in the dissipation oc- 
curs at a temperature-dependent field Hj(T) (Fig. 2), 
which we associate to Josephson vortices entering the 
sample. Interestingly, above Hj(T), the microwave ab- 
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FIG. 2. Dissipation as a function of the applied field at 3 
temperatures, for 2 orientations (0° and 2°) of the applied 
field with respect to the afo-plane. The onset of dissipation, 
indicated by the arrow, occurs at the penetration field Hj(T). 



sorption behaves linearly with field, with a very good 
accuracy, from typically 10 Oe up to 25 Oe. This ap- 
pears consistent with a flux-flow mechanism driven by 
c-axis currents, where the flux-flow resistivity is linear 
with applied field. We therefore identify Hj(T) to the 
first penetration field of Josephson vortices. In this work, 
unlike in Ref.EZI, we have averaged the data over the field 
orientations from 0° to 3° , in an attempt tp-improve the 
accuracy of the determination. As in Ref£3, we choose 
to define Hj(T) as the field value where the microwave 
absorption exceeds the experimental accuracy (2 • 10~ 4 ). 
The field thus determined is plotted in Fig. 3. The error 
bars take into account both-thc noise and the estimated 
drift of the signal with timea. 



IV. DISCUSSION 

Whether Hj(T) may be identified to the ther- 
modynamic— lower critical field H c \ was previously 
discussecOlEI. For Josephson vortices H c i writes^: 



H cl {T) 



4ir\ ab (T)\ c (T) 



[lnA a6 (T)/d+1.12] (1) 



In our early workE3, we had not yet studied the temper- 
ature dependence of Hj(T) and we had not observed the 
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FIG. 3. Plot of Hj(T) (full circles). Up (down) triangles 
display estimates of H c i(T) using A c (0) = 10 (jm (40 /im). 
The temperature variations AX a b(T) and AX C (T) are taken 
from our present work (Figs. 1,5). 



irreversible behavior of the dissipation upon flux entry 
and flux exit. We had therefore not considered the possi- 
bility of a surface barrier. However, in order to reconcile 
the magnitude of Hj(T << T c ) with the thermodynamic 
lower critical field, we were compelled to take the lowest 
possible values for A a &(0) and_A c (0). 

We proposed next in Ref.E3 to take more acceptable 
lower bounds for A a fc(0) and A c (0), together with the ex- 
perimentally determined temperature variations in order 
to obtain an upper bound for j? el (T.) .. Here, we take 
2100 A as, r-a, iawer bound for A a6 (0)ElH, and 10 /im 
for A c (0)E3I23e3. We use the temperature dependence 
for AA a fc(T) (partly shown in the inset, of Fig. 1) and 
AA C (T) measured in our previous workn3 (see Fig. 5 be- 
low) . The corresponding H c i (T) is plotted in Fig. 3 using 
the above mentioned values. We havei a^so displayed in 
Fig. 3 H cl (T) if taking A c (0) = 40 /imBE3. It is clearly 
seen that neither the absolute value (too small compared 
to the experimental data) nor the temperature depen- 
dence (quasi-linear) agrees with the Hj{T) data. Since 
the actual penetration field is larger than the thermo- 
dynamic H c i(T), it is therefore quite natural to assume 
that a Bean-Livingston surface barrier impedes field pen- 
etration, and yields a larger entry field H^ D (T). 

In anisotropic superconductors, in the quasi-2D 
regime, e.g. when the transverse coherence length £ c be- 
comes smaller than the interlayer distance d, H^ D (T) was 
shown ta be related only to the c-axis penetration length 
througljj : 



H 2 e D {T) = 



An\ c {T)d 



(2) 



In Bi-2212, the quasi-2D regime holds up to temper- 
atures very close to T c , so that this last expression for 
H^ D (T) is valid in our measuring temperature range. A 
surface barrier might thus account for the observed value 
of the penetration field. Also, since H^ D (T) grows as 
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1/A C (T) (instead of l/A a &A c (T)), it is expected that the 
temperature dependence could show a better agreement. 
The existence of a surface barrier is further suggested by 
the hysteretic behavior of dissipation, shown in Fig. 4, 
at T — 65 K (the behavior is similar at other temper- 
atures). When the field is swept down, vortices do not 
exit in a reversible way. However all vortices have left 
the sample as can be inferred from the recovery of the 
same dissipation as in zero initial field, when the field 
is back to zero value. When the field is swept up again, 
the absorption displays precisely the same behavior as af- 
ter the zero field cooled procedure. Bulk pinning would 
induce flux trapping at zero field, hence some residual 
dissipation. Our observations are similar to magnetiza- 
tion measurements where the irreversibility, assigned to 
a surface barrier, is characterised by zero magnetization 
upon decreasing field. Such apbehavior, first observed 
in a field parallel to the c-axisE-3, was also reported for 
Josephson vortices in Bi-2212 in a field oriented nearly 
parallel to the ab planets—Surface barriers may also lead 
to time dependent effectscEl. Indeed, it was argued from 
magnetization data taken at various sweep rates that the 
penetration field in the parallel configuration depends on 
the field sweep rate, and eventually achieves the thermer 
dynamic first critical field value for the slowest ratesc3. 
Our sweeping rate is of the order of 0.1 Oe/s, compara- 
ble to the_range where the largest penetration fields were 
observecO. We did not change the sweeping rate hence 
we cannot confirm this claim. We point out however 
that the penetration fields observed ir£3 are significantly 
smaller (roughly a factor of 3) than ours. Compared to 
the fastest rate, the decrease of the penetration field asso- 
ciated to the slowest rate is only 1 Oe. Such small values 
can obviously be more easily reconciled with H c \ (T) than 
ours. 
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FIG. 4. Plot of x"(H) as a function of the direction of the 
field sweep. Full circles, full squares and open squares refer 
to sweeping up the field starting from a ZFC state, sweep- 
ing down the field to zero, and sweeping up the field again 
respectively. 

It is worth noting that all these remarks do not modify 



the surface barrier interpretation: they only put a time 
scale for its observation. 

Relying on the results described above, we derive from 
the Hj(T) data an effective penetration depth Aj(T) us- 
ing Eq. (2). The data are shown in Fig. 5. We then try 
to determine A c (0) so as to fit Aj(T) using the measured 
AA C (T). We find that both sets of data, namely Aj(T) 
and AA C (T) cannot be reconciled for any value we may 
assume for A c (0). Therefore, the interpretation cannot 
be so simple. 
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FIG. 5. Plot of the temperature variation AA C (T) (full cir- 
cles, right scale) and of the effective length Xj(T) which is 
associated to a surface barrier for the penetration of Joseph- 
son vortices (full squares, left scale), using Eq. (19) or (21). 
Open circles display the best fit using A c (0) = 35 /im and a 
scaling factor (3 = 6. Open symbols show the best fits using 
A c (0) = 20 /im (down triangles) and A c (0) = 100 /jm (up 
triangles). 



V. ROLE OF SURFACE IRREGULARITIES 

A. Experimental checks 

Actually, a surface barrier is only effective if the surface 
is smooth on a typical length scale which is the penetra- 
tion depth. In our field geometry, defects either located 
on the top and bottom, e.g. a6-planes, or the edges may 
destroy the surface barrier. In the former case, the rele- 
vant length scale is A a &(T), in the latter case, A C (T). 

In order to distinguish between these two possibilities, 
we have carried out several checks. The samples dis- 
cussed in this paper were first measured without any spe- 
cial preparation except for their initial shaping in platelet 
and cleaving in order to work on a well defined single 
crystal and mirror-like surfaces. We noticed that cleaved 
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surfaces often exhibit a few visible steps and sparse voids. 
After the first measurement, sample D was placed on 
the stage of an STM and the tip was used in order to 
cut four grooves parallel to the small side of the crystal, 
4000 A deep and 100 /im apart. Then the sample was 
measured again (H parallel to the grooves). No signifi- 
cant change in the onset field of the microwave absorp- 
tion was observed. In a second step, we took another 
sample yielding a similar penetration field, and cleaved 
it. We obtained fresh surfaces with one or two isolated 
steps which could be seen under a binocular. This sample 
was measured immediately after cleaving, and again, no 
significant change was observed. It seems therefore that 
either defects within the afr-planes do not play any role 
in order to reduce a surface barrier or even a single step 
is immediately effective to destroy the surface barrier. 

One should also consider penetration through the 
edges. Indeed, the edges are fairly difficult to control. 
We did check indirectly, in the surface impedance and ac- 
susceptibility experiments, whether they play any role. 
In order to measure AA C (T), the rf magnetic field ap- 
plied parallel to the plane is also parallel to one edge 
of the crystal. If the sample is rotated by 90° along 
its c-axis, the edges where c-axis currents flow are in- 
terchanged. It is then clear that if there exists a large 
defect, e.g. a slit or groove deep in one edge and not 
in the other, this defect changes significantly the c-axis 
microwave current pattern in one position and much less 
in the other. Therefore, the two configurations should 
yield a different AA C (T) result. In one particular sam- 
ple out of three, this was indeed the case, suggesting the 
presence of a defect lying in one edge and showing that 
the measured AA C (T) cannot be intrinsic for this partic- 
ular sample. The data used in this paper and shown in 
Fig. 5 are not biased by such large edge defects, e.g. the 
AA C (T) data displayed in Fig. 5 are unchanged within 
the accuracy of the measurement upon this rotation. 

We have now to examine quantitatively to which extent 
defects located on the top or bottom surfaces, or in the 
edges alter the penetration field. 



factor of two near the defect. Then, the main change of 
the entrance field is essentially related with the increase 
of the screening current density. 

We consider the case where the scratch is in the form 
of a groove on the superconductor surface, and the mag- 
netic field is parallel to it. Let the z-axis be perpendicular 
to the superconductor surface. The magnetic field is par- 
allel to this surface along the y-axis, and we choose the 
axis of the groove on the same direction. The depth of 
the scratch is denoted as b and its width 2a (see Fig. 6). 
For convenience, the semi-axis z > 0, is chosen inside 
the superconducting material, so in Fig. 6 the scratch is 
presented on the bottom surface of the superconductor. 
Both a and b are supposed to be much smaller than A, 
the London penetration depth, so screening can be ig- 
nored and the two-dimensional London equation reduces 
to Poisson's equation. Then the lines of current corre- 
spond to the equipotentials, and a dielectric defect in a 
superconductor corresponds to a metallic embedding in 
electrostatics^. This analogy reduces our problem to the 
calculation of the electric field distribution near a metallic 
electrode having the special form (Fig. 6) while the field 
becomes uniform for z-~-> oo. As is known from electro- 
statics (see, e.g., Ref.c3), the solution is provided by a 
conformal transformation of the w-plane, corresponding 
to a flat surface, to the £-plane, the plane of an orthogo- 
nal cut of the scratch. In the w-plane the attraction en- 
ergy of the plane on a vortex located at the point w can 
be easily computed, for example by the image methodo 
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B. Theoretical calculations: formalism 

As usual the entrance field is deduced from the bal- 
ance between the vortex attraction to the surface and 
the pushing force exerted by the screening-current at the 
minimum distance £ (the vortex core size)Bc3. The pres- 
ence of the surface irreqularities can strongly influence 
the screening current distribution. In particular, near a 
scratch the current density can be many times larger than 
near the flat surface. This may substantially increase 
the force pushing vortices inside the superconductor and 
then decrease the surface barrier and the entrance field. 
The vortex attraction to the surface does not change es- 
sentially near a scratch, as it has been demonstated in 
Rcf.c3. The force of attraction can decrease by at most a 
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FIG. 6. a) The defect MABCDN in the form of the groove 
at the bottom surface of superconductor. The plane of the 
figure corresponds to the plane £ = x + iz. b) The plane 
w — u + iv where the straight line M'A'B'C'D'N' is mapped 
to the surface line MABCDN. 
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Besides, a uniform current density j(w) = jo in the w- 
plane, can be deduced from the simple complex potential 
ip(w) = jotv. Then in the £-plane, the complex current 
density j (£) can be obtainded from the complex potential 
■0(C) = jo w(() by 
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(4) 



the entrance field. In fact this situation where z ~ £ ^> a 
is not realistic for isotropic superconductors, but it will 
be useful for the description of the anis otrop ic ones 



For the groove-like defect, equations ( A.14 ) and ( A.16 ) 
allow to derive the physical quantities in the vicinity of 
the point C. Let P be a point such that 



Cp = a + pe , p < a. 



(9) 



where jo is the current density far away from the defect, 
i.e., the screening current near the surface jo — cH/AnX. 
To calculate both the attraction energy and the current 
density in the £-plane, we need to inverse the conformal 
transformation. In general, this cannot be done analyt- 
ically. However, for situations of practical interest, we 
may use approximations that allow us to obtain an ana- 
lytic solution. 



C. Isotropic case 



The values of 8 are limited by the groove and the core 
of the vortex 



arcsin(£/p) < 9 < ir — arcsin(£/p). 



(10) 



The attraction energy on a vortex at the point P is 



E att (p, 0) 
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In the Appendix, we have demonstrated that according 
the values of a, b, £, there are three different regimes: 

(i) a < b — |C| slit-like defect, 

(ii) |C| <C a <C b groove- like defect, 

(iii) |C| <C b <C a step-like defect. 

Let us begin by th e slit -like defect. In this case, we use 
equations ( |A9| ) and ( A12 ) to derive the vortex attraction 
energy at the distance z from the slit 
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and the strength of the attraction force is 



fatt(z) 



2 V4ttA / z 



(5) 



(0) 



that is half of the force for a plane surface. Similarly the 
current density is 



]{Z) = T2 30 



(7) 



These two results were obtained earlier0 for the field 
and current distribution near the angle 2n (the cut at the 
superconductor surface). 

The balance between the vortex attraction to the sur- 
face and the pushing force exerced by the screening aic. 
rent at the minimum distance £ (the vortex core size)l3'E^I 
gives the entrance field near the defect 



H P . 
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(8) 



and the strength of the attraction force reads 

2 



fatt(p,0) 



= 2 (4o_\ 1_ 
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(12) 



Its maximun is obtained for 9 = or 7r/2, this strengh 
is reduced by a factor 4/(3V3) — 0.77 by comparing to 
a flat surface. The calculation of the current density at 
the point P gives 
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As usual by setting a vortex at a distance £ near the 
defect, we obtain for the entrance field 
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(14) 



where [16/(90F)| 1 / 3 ~ 1. 

Fin ally, for a step a 3> b, by using equations (A14) and 
( A17) we can derive the vortex attraction energy and the 
current density at the point P in the vicinity of the point 
C. The vortex attraction force is still given by Eq. jl^ ) 
while the current distribution near the corner becomes 



1/3 



b 

JO \ ~ 

P 



1/3 



-i0/3 



The corresponding entrance field is 
1/3 / c \ 1/3 
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(15) 



(16) 



where H e is the entrance field for the flat surface H e ~ 
4>o/iirXt; ~ H c (thermodynamic critical field). The cur- 
rent concentration effect near the slit essentially reduces 



where (27T/9) 1 / 3 ~ 0.89. 
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D. Anisotropic case 

Now we consider the case of anisotropic superconduc- 
tors, keeping in mind layered high-T c materials. As 
usual, let the z-axis (or c-axis) be perpendicular to the 
superconducting layers. We shall consider two cases : ei- 
ther the groove is on the bottom surface of the crystal, 
or it is on the side surface (edge) of the layered material. 
In both cases we choose the axis of the groove parallel 
to the layers along the y-axis, and the magnetic field in 
the same direction : h=h(x,z)e y . For such a geometry 
we may write the London free energy of the anisotropic 
superconductor as 



Near the step a' 3> b 3> d, the entrance field becomes 



F 



1 

8^ 



h 2 + 



dh 
dx 



+ A 



dh 



dV, (17) 



where A c is the London penetration depth when the 
screening current is flowing along the z-axis (c-axis) and 
X a b when the current is in (x,y) plane. For a high-T c 
superconductor, we have A c ^> A a &- 

For a very anisotropic superconductor, in the quasi 
2D regime we have £ c < d, where d is the interlayer 
distance. In such a case, we need t©, use d as the size, 
in the z-direction, of the vortex coreQ in calculating the 
entrance field. For the flat surface the entrance field be- 
comes H 2D ~ o /[47rdA c ]. 

By making a scaling transformation, we, introduce a 
new coordinate : X — (X a b/X c )x -C xlH. Then the 
London free energy ([l?]) takes the same form as for the 
isotropic superconductor with the London penetration 
depth A a fc and we can use the results of the previous 
section. 

Let us consider the case when the groove is on the 
bottom surface of the crystal. Under the scaling trans- 
formation, the width of the groove changes 



A, 



■ a <C a. 



(18) 



Then, the entrance field will be given by the corre- 
sponding formulas for the isotropic case with the replace- 
ment of a by a' and £ by £ c , the correlation length along 
the z-axis (or by d when £ c < d). 

For d ~ b 3> a! , the groove may be considered simply 
as a thin cut at the surface and by using Eq. (ph we derive 



tt2D 



H. 



2D 



1/2 



(19) 



Note that due to the large anisotropy of some high-T c 
materials, this situation could be realized in practice de- 
spite the fact that d is of the order of only 10 A . 

In the opposite cased <C a' <C b, Eq. ( [l4| ) gives the 
entrance field near the groove 



tt2D 



H 



2D 



1/3 



aXg, 
~bX c 



1/6 



(20) 



H 2D ~ TT 2D 



1/3 



(21) 



Finally, if the groove-like scratch is on the side surface 
of the layered material, its effective depth b' after the 
scaling transformation becomes much smaller 



b' = ^ b « b. 
X c 



(22) 



Then for b' <C a such a scratch has practically no effect 
on the vortex entrance. In the opposite case b' 3> a > d, 
by using Eq. (fTi]), we obtain for the entrance field : 

In fact the lateral defect must be rather deep: b ^S> 
aX c /Xab > dX c /X a i> to strongly reduce the entrance field 
value. We may deduce that the parallel entrance field 
depends strongly on the surface defects in layered super- 
conductors. The current concentration near the defect 
edges may greatly reduce the entrance field in compari- 
son to its theoretical value H 2D ~ o /47rdA c . 



E. Comparison with experimental data 

We have therefore attempted to fit our data derived 
from Hj{T) with Eqs. (19), (21) or (23), using two ad- 
justable parameters: a scaling factor /3 associated with 
the defect geometry which equals to (b/d) 1 ^ 2 in Eq. (19), 
(b/d) 1 / 3 in Eq. (21), (d 2 a/b 3 ) 1 / 6 in Eq. (23), and the ab- 
solute value of A c (0). We show the results in Fig. 5, only 
for the case described by Eqs. (19) and (21) (defect in 
the a&-plane) , where we have determined the scaling fac- 
tor and A c (0) which allow to adjust Xj(T). We find a 
best fit for A c (0) = 35 ^m and a scaling factor /3 = 6. 
Assuming a thin groove, this yields b ~ 500 A which is 
reasonable. We also show in Fig. 5 smaller and larger 
values for A c (0). They allow us to set the uncertainty 
about our determination of the penetration depth. 

As for Eq. (23) (edge defect), we have taken A a fc(0) = 
2600 A . We can also account for the data but only in a 
very restricted, nevertheless acceptable, range of param- 
eters. The depth of the edge slit should be of the order 
of 10 /im, which is still small with respect to A c (0). The 
key result in this latter case is that it yields the same 
absolute value for A c (0). 

In conclusion, the set of experiments that we have per- 
formed suggest very strongly a surface barrier which im- 
pedes field penetration, nevertheless partially destroyed 
according to the calculations developed in the framework 
of this work. Although we cannot ascertain which spe- 
cific defects reduce the efficiency of the surface barrier, 
we obtain a fairly good estimate of the c-axis penetration 
depth. 
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APPENDIX: CONFORMAL TRANSFORMATION 

The Schwarz-Christoffel conformal transformation of 
the w plane (w = u + iv) to the C plane (£ = x + iz) 
which maps the straight line M'A'B'C'plN' (Fig. 6b) to 
the surface line MABCDN (Fig. 6a) isN 



CM = A 



u/e 



l-t 2 
1 - k 2 t 2 



dt. 



(Al) 



where the two parameters k and £ are related to the di- 
mensions a and b of the groove, and the constant A is 
simply 



A' 1 =a 



1-t 2 
1 - k 2 t 2 



dt. 



(A2) 



The integrals of the two previous equations can be ex- 
pressed in terms of incomplete and complete elliptic ir 
grals E(z, k), F(z, k), E(k) = E(l, k), K(k) = F(l, fc)E 
For this we define two G functions, one incomplete and 
one complete as : 



G(z,k) = E(z,k)-k' 2 F(z,k). 
G{k) = G(l, k)=E (fc) - k' 2 K (k) . 



(A3) 



where k' = \/l — k 2 . Then the conformal transformation 
reads 



CM = « 



G(w/e, k) 

G{k) 



(A4) 



The dimensionless parameter k 6 [0,1] is determined 
by the following equation, 



G(k) 
G(k') 



(A5) 



The limits k — > and k — > 1 correspond to a/b — > 
and a/b — > oo respectively. The last parameter £, the 
dimension of which is a length, is determined by requiring 
that at a large distance from the defect, the two variables 
C and w become equal. Using the asymptotic form of 
G(z, fc) for large z, 



£ = a 



G(k) 



(A7) 



When k — > 0, i.e., k' — > 1, the following asymptotic 
forms of the elliptic integrals 

7T , 



0, G(fc) ~ -fc 2 , G(fc') ~ 1, 



(A8) 



may be used to determine the parameters k and £ in the 
limits a/b — ► or a/b — > oo 



a/b^0, k: 
a/b — > oo, k ~ 1 



;t V b 
2 b 
7r a ' 



'a b . 



(A9) 
(A10) 



Firstly, we suppose that the groove is very narrow a <C 
6, and we consider the region where C ~ w ~ b. Then we 
have \w\ 2> £ and the elliptic integrals can be approximcd 

as 



k -> 0, 1 < |z| < 1/fc G (z, fc) ~ 2 fc ^ AU ) 

The conformal transformation becomes simpler and it 
can be inverted 



CM 



l w 
2T 



(A12) 



Secondly we consider the vicinity of the point G in the 
C-plane (|C — a\ a) and of C in the w-plane (\w — l\ <C 
£), In this case we have \w/£— 1| <C 1 and the behaviour 
of the elliptic integrals is 



\ v \<l^G(l + r),k)~G{k)-i 



2%/2fc 2 
3fc' 



3/2 



(A13) 



As previously the conformal transformation can be eas- 
ily inverted 



CM - a - -?= 



3/2 



fa \ (p(k) 
w(Q ~ £ = e"/ 3 (^(fc)a 1/3 (C - a) 2/3 , 
where the function ip(k) is defined as 



<p(k) 



9fc' ; 



1/3 



kG (fc). 

The asymptotic forms of this function read 



a/b^O, ^^(ffl 



1/3 



a/6 oo, p(fc) ~ ( — 



9\ 1/3 /6 



(A14) 



(A15) 



(A16) 



(A17) 



Note that in this last limit a ^> b, we retrieve the case 
of a single step defect. The second step at the point G is 
not involved. 



\z\ > 1/fc => G(z, fc) ~ kz, 



(A6) 
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